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Abstract 

Using a Magnetic Resonace Force Microscope, we have performed ferromagnetic resonance 
(FMR) spectroscopy on parametric magnons created by 4-wave process. This is achieved by mea- 
suring the differential response to a small source modulation superimposed to a constant excitation 
power that drives the dynamics in the saturation regime of the transverse component. By sweeping 
the applied field, we observe abrupt readjustement of the total number of magnons each time the 
excitation coincides with a parametric mode. This gives rise to ultra-narrow peaks whose linewith 
is lower than 5 10~ 6 of the applied field. 



The detailed understanding of the non-linear (NL) regime of the magnetization dynamics 
is important both from a fundamental point of view-i but also for applications in spintronic 
devices^. Interest resides in the exact nature of the parametric modes that are excited above 
the supercriciality threshold. Recent experiments performed on a yttrium iron garnet (YIG) 
film have shown that these parametric magnons can form a Bose-Einstein condensate^ under 
high power pumping. It was also demonstrated that their energy decay rate to the lattice 
is substantially diminished compared to the long wavelength modes^, which are usually 
studied by FMR In this paper, it will be shown that it is also possible to measure the FMR 
spectrum of these parameric magnons despite the fact that these high /c-vector spin-waves 
do not overlap with the homogeneous microwave field. This is achieved by saturating the 
dynamics of the mode that couples to the microwave. Any additional power that is injected 
in the spin system flows directly towards the excitation of these parametric magnons. Since 
each parametric mode has a different feedback influence on the saturation, it is possible, by 
sweeping an external parameter, to detect the point where the transfer from one parametric 
mode to another takes place. But the induced effects are small and detecting them requires 
both a sensitive and precise measurement setup. 

Our detection of the dynamical response uses a mechanical-FMR setup working at room 
temperature^. This scheme is inspired by local probe techniques. The static component 
of the disk magnetization, M z = M s cos8, is coupled through the dipolar interaction to 
a magnetic probe attached at the end of a soft cantilever, whose deflection is monitored. 
Exciting the sample at a fixed frequency, spectroscopy is achieved by recording the cantilever 
motion as a function of the perpendicular dc applied field, H ext , which is produced by a NMR 
electromagnet. The variation of force on the cantilever is proportional to the total number 
of magnons excited at resonance, N t , with N t = (M s — M z ) / (jh). Further details about the 
mechanical-FMR experiment can be found in Ref.-. 

Driving the dynamics of ferromagnets in the NL regime requires to cant the magneti- 
zation, M s , away from its equilibrium axis by an angle 9 exceeding a couple of degrees. 
Only powerful HF source can produce such large excitation amplitude^. Improvement of 
the efficiency can be gained with microstrip cavities, where the HF energy is concentrated 
in smaller volume. A schematic of our setup is shown in FigfjTa). It comprises a 0.5mm 
wide Au stripeline fabricated by U.V. lithography and deposited on a 0.5mm thick alumina 
substrate whose bottom layer is a conducting ground plane. An impedance matched cavity 
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(half-wavelength) is created by etching a 32/im gap across the stripe. The HF source is 
tuned at the cavity frequency (u s /2n = 10.47GHz) set by the length (5mm) of the isolated 
segment. The sample is a t = 4.75/im thick YIG single crystal. The YIG film is ion milled 
into a disk of diameter <fi = 160/im and placed at the center of the cavity. The homogeneous 
external static field H ext , well above the saturation field, is applied parallel to the disk axis. 

This work concentrates on the 4-magnons coupling term in the equation of motion of 
the magnetization. It becomes the dominant term once the premature saturation regime of 
the microwave susceptibility is achieved^. From previous experiments on the same sample^, 
we have established that the saturation threshold occurs at excitation power above -35dBm 
for our setup. Figfjjb) shows the static deflection of the cantilever as a function of H ext 
when P max = —10.8 dBm. The trace corresponds to the lineshape of the uniform mode^ 
in the foldover regime. Hysteretic behaviors are a standard signature of NL effects, where 
the eigen-frequency of the resonance depends on the amplitude of the excitation^. Foldover 
effects in magnetic materials have been explained by Anderson and Suhl in 1955^. The 
features of interest are the tiny steps observed in the wings of the resonance (see insert of 
Figfjjb)). These jumps are not random. They are reproducible and occur regularly in field. 

A more detailed picture can be obtained by measuring the differential response. We 
detect here the response to a source modulated excitation around the saturation regime. 
The power level follows the time dependence P(t) = P max {l + § (sin(u; c t) — 1)}, where the 
modulation frequency is set at the resonance frequency of the cantilever. This frequency 
is much lower than all the damping rates of the spin system. In this case, the differential 
part of the mechanical signal is amplified by Q = 4500, the quality factor of the mechanical 
resonator. Figfjjc) is the pattern recorded by a lock-in when the modulation amplitude is 
0.5% (e = 5 10 -3 ) of the maximum power, P max = —17 dBm. It corresponds to a source 
modulation of less than 50nW in amplitude. The spectrum shown here is a zoom on a 2G 
sweep of ifext in the reversible region of the static signal. We observe about ten ultra-narrow 
lines (width 30mG) regularly spaced every 60mG. Such sharp peaks resemble the parametric 
excitations observed by Jantz and Schneider in 1975 in the subsidiary absorption of YIG 
films^ 3 - (3-wave process). To the best of our knowledge, this is, however, the first time that 
they are observed at resonance. 

The striking characteristics of these peaks is their narrowness. This 30mG broadening 
is much smaller than the FMR linewidth AH = 1.2G measured in the linear regime on 
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the same sample^. Further insight can be obtained by looking at the shape of the signal 
(see FigfS^) for different modulation depth, e, at constant maximum power P max . Figfjb 
shows that the broadening B of the parametric peaks is clearly proportional to e. Foldover 
effects establish a coupling between the HF power and the resonance frequency. Thus source 
modulation in the NL regime produces also a frequency modulation. But in FMR, there 
is no difference between frequency modulation and a modulation of the applied field, since 
both are coupled by the gyromagnetic ratio. Fig{Ib can thus been seen as the differential 
characteristics of Figfjb, where the depth of modulation is proportional to e. At very low 
e eventually the intrinsic broadening should be seen. Although a NMR electromagnet with 
a 10 -6 precision in if ex t has been used, the fact that B extrapolates linearly down to the 
smallest value of e reached by our setup, suggests that the intrinsic part of the broadening 
has not been measured here. 

For completeness, we show in FigfS]how the instability region evolves when the dc bias 
P max is increased at constant e. We observe an increase of the applied field window, where 
parametric peaks occur. The size of the window follows monotically the amplitude of the 
dc power i.e how far we depart from the supercriticality threshold. 

In the following, we propose an analytical framework to account for the main features 
observed experimentally. Saturation phenomena can be caused by a multitude of nonlinear 
processes: nonlinear dissipation, nonlinear phase mechanism, nonlinear energy feedback 
effect. In many cases several of these mechanisms, as well as nonlinear frequency shift, need 
to be accounted simultaneously for the correct description of the phenomenon. The most 
probable limiting mechanism in the case of perpendicular pumping is the nonlinear energy 
feedback. This mechanism is described by a 4-magnons interaction (Cfc^o^-fc + c - c ) i n the 
Hamiltonian of the system. It represents second-order parametric excitation of the pair of 
plane spin waves (bk, 6-fc) by the uniform precession (6 )- Suhl showed in 1957, the equations 
of motion for the amplitudes bk of the plane spin waves have the form: 

b k = -iuj k b k - i£kbo 2 b- k * - f] k b k , (1) 

where rjk is its energy decay to the lattice (the index must allow for fc-dependent relaxation 
rates^i£) . The equation for the amplitude bo of the uniform precession is given by 

b = -iu b - % ikbkb-kbl - Vob + -fhe'^ 1 . (2) 
k 
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where the last term is the excitation field h = \f~P '. The term Y2k ■ • • describes nonlinear 
feedback effect of all the parametric waves on the uniform precession. This term gives 
effective nonlinear dissipation for 6 - m t ne following, we will be looking at harmonics 
solution where all the spin- waves precess synchronously: b^ = —iuj s bk- 

As one can see from Eq. (TjQ), an instability (exponential growth of the population) occurs 
if the effective damping of at least one parametric mode becomes negative: 

,4 . (gfc ~ u s ) 2 + rg 

l & 0| > 72 • K 6 ) 

Growth of bk will create effective damping for the unifrom mode 60 (through the Ylk • • ■ 
term in Eq. (j2j)), which will reduce 60 to its threshold value: 



\h 12 AT - ■ - Us) 2 + r£ f .s 

\b Q \ = N = mm , (4) 

k gt 

Determination of the mode that will grow unstable depends on the coupling which is 
maximum ^k/l — 27iM s (^i 900G) for k propagating parallel to the magnetization direction^. 
As shown on the magnon manifold drawn in FigJU these longitudinal SW have a wave- 



vector A; max = y/N ± 4irMj/(2A) w 6.3 10 4 cm , where A is the exchange constant and iVj_ 
is the transverse depolarization factor. For our disk, it corresponds to a standing spin-wave 
confined across the thickness with about 5 nodes along it. In this model, the threshold power 
is given by the analytical formula, h 2 = r]kAH 2 /C,k and the calculated value h c = 5 mOe 
agrees with the the onset of saturation measured experimentally 4 -. 

Sweeping the bias magnetic field if ox t (or u s ) changes the nature of the parametric mode 
that grows unstable, as the minumum in Eq. (j4j) occurs for different value of k. Finite 
size-effects introduce a discretisation in fc-space of all the SW. The frequency separation Au 
between two nearby normal mode is about Au; = 4^AkAk/M s , where Ak = n / cf) corresponds 
to the quantification of the fc-vector along the largest sample dimension. The parametric 
modes have thus a principal wavector k max \\M s , with a small additional component A A; 
in the disk plane. We find numerically Acj/7 = O.120e, a predicted separation twice as 
large as what is observed experimentally in FigJT](c). The value of /c max however depends 
strongly on the angle of propagation of the degenerate magnons. A deviation of about 7° of 
their propagation direction compare to the normal of the disk will be enough to explain the 
discrepancy. 

The amplitude of the steps can also be estimated analytically with an approximated 
model. Assuming that only one pair of parametrically coupled waves dominates all the 



others, then its amplitude can be infered from Eq. (jSJ). If the mode b k is excited, then its 
amplitude will grow to the level: 



\ (jh) 2 ^iVo - {(u k - u s )r) + (u - u s )Vk} + (^o - u s ) (u k - u s ) - r) r) k 
\b k \ 2 = \bA 2 = N k = -i = 

N e k 

(5) 

that provide the effective damping to fix iVo at the level given by Eq. (J4]). The number of 
magnons in the majority mode N k changes abruptly with the change of the excited mode due 
to the sweep of the magnetic field. In the differential graphs this will lead to a singularity at 
the points, where the transfer from one excited mode to another takes place. In reality, due 
to the presence of the thermal noise and other nonlinear processes (that are ignored here) 
one will observe a sharp (intrinsically broaden) peak at the points where the transfer from 
one excited mode to another occurs. 

To see this more clearly, let us assume that coupling coefficients of all modes are the 
same, = £. Then the changes of the parametrically excited modes (say, modes 1 and 2) 
will be at such a magnetic field, that (oui — u s ) = — (a^ — uj s ) = Au/2 (see Eq. (j4j). 

Assuming that \uq — uj s \ <C 77^ and Au; -C rjk we can write iVo ~ ?7fc/£ and, approximately, 

^(c-Df ± (C- Dfo- (6) 

where ( = h/h CT is the supercriticality parameter of the microwave magnetic field. 
Thus, the total number of spin wave modes differ by 

which will give small "steps" (of the order of 1% of N t ) in the behavior of M z (H cxt ) at 
different "integral" characteristics and ultranarrow peaks at "differential" characteristics 
of M z vs. H cxt . These analytical findings are in good agreement witht the experimental 
observed in FigJTJ 

In conclusion, we have achieved FMR spectroscopy of parametric magnons. The ob- 
served linewidth of 30 mG corresponds to a lifetime of the order of 2/is, which should be 
compared with the spin-lattice relaxation rate T\ = 106ns usually observed in the lineare 
regime^. Such long lifetime starts to become comparable to the one found in nuclear mag- 
netic resonance using shimmed magnets. This result opens up the possibility of making 
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ultra-stables Yig Tuned Oscillators or ultra-sensitive absolute magnetic field detectors using 
these parametric resonances. 
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FIG. 1: (a) Schematic of the microstrip circuit (1). The FMR spectrum of a micron-size YIG 
disk (2) is detected by Magnetic Resonance Force Microscopy (3). On the right is shown P(t), 
the waveform of the power injected in the microwave circuit (1). The amplitude of the sinusoidal 
modulation is eP max /2. (b) Static deflection of the cantilever as a function of H e xt- The trace is 
the lineshape of the uniform mode in the foldover regime, (c) Vibration amplitude of the cantilever 
measured by a lock-in synchronous with source modulation. The differential characteristics reveals 
the parametric excitation spectrum. 
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FIG. 2: (a) Differential spectrum as a function of the modulation depth e (constant P I 
— 17dBm). B is the broadening of the parametric peaks, (b) B as a function of e. 
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FIG. 3: Differential spectrum as a function of P max (constant e = 0.1). 
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FIG. 4: (a) Magnon manifold of our perpendicularly magnetized YIG disk with a schematic of 
the mode discretisation induced by the confinement in a micron-size sample. Four-waves process 
couples the uniform mode (red) to the parametric modes at 9^ = (green) (b) Schematic of the 
distorsion induced by the 4-waves NL process on the static spectroscopy of M z = 7^X/fe-^fe vs - 
u s . The staircase behavior reveals the parametric spectrum in the differential characteristics. 
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